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Abstract 



This paper is the first in a project concerned with the quantum- 
mechanical decay of a Schwarzschild-lilce black hole, formed by gravita- 
^ , tional collapse, into almost-flat space-time and weak radiation at a very 

QQ ' late time. The approach taken here is, in fact, applicable to a much 

I wider class of quantum calculations than those concerning gravitational 

, collapse, including quantum properties of a variety of cosmological exam- 

' pies involving weak anisotropic perturbations. In this work, we are con- 

cerned with evaluating quantum amplitudes (not just probabilities) for 
. transitions from initial to final states. The present quantum description 

shows that no information is lost in collapse to a black hole. In a general 
^ , asymptotically-flat context (not necessarily involving local gravitational 

O^' collapse to a black hole), one may specify a quantum amplitude by posing 

boundary data on (say) an initial space-like hypersurface E/ and a final 
{3/), space- like hypersurface Ef , together with the Lorentzian proper-time in- 

terval T which separates them, as measured near spatial infinity. Suppose, 
for simplicity, that the Lagrangian contains Einstein gravity with only a 
■ minimally-coupled massless scalar field , describing the matter present. 

' Then the boundary data can (for example) be chosen to be hij and (f} 

on the two hypersurfaces, where hij = 1,2,3) is the intrinsic spatial 
3-metric. The classical boundary-value problem, corresponding to the 
calculation of this quantum amplitude, is badly posed, being a boundary- 
value problem for a wave-like (hyperbolic) set of equations. Following 
Feynman's +ie prescription, one makes the problem well-posed by rotat- 
ing the asymptotic time-interval T into the complex: T \T\ exp{—i9), 
with < 9 < tt/2 . After calculating the amplitude for 6 > , one then 
takes the 'Lorentzian limit' 9^0+. All calculations in this work are 
based on this procedure. For example, in the following Paper II, this will 
be used to calculate amplitudes for given final weak configurations of the 
scalar field, representing scalar radiation on the final hypersurface Ef at 
a late time T . 
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1 Introduction 



This paper is the first in a project concerned with the calculation of quantum 
amplitudes (not just probabilities) associated with quantum fields (including 
gravity), in the case that strong gravitational fields may be present; this may be 
compared with the effective-action approach of [1,2]. The most obvious example 
- the original motivation for this work concerns quantum radiation associated 
with gravitational collapse to a black hole [3-11]. But the framework adopted 
here is more general, and certainly does not depend on whether there is a clas- 
sical Lorentzian-signature collapse to a black hole. It includes the case of local 
collapse which is not sufficient to lead to (Lorentzian) curvature singularities, 
and also quantum processes in cosmology, where, for example, anisotropics in 
the Cosmic Microwave Background Radiation can be computed, and depend 
crucially on the underlying Lagrangian for gravity and matter [12]. 

To exemplify the underlying ideas, let us consider the case of local collapse 
(whether or not to a black hole). Thus, the gravitational field is taken to be 
asymptotically flat. Suppose, as in this paper, that we consider Einstein grav- 
ity coupled minimally to a massless scalar field <j> . In classical gravitation, we 
are used to describing this by means of a Cauchy problem, giving evolution to 
the future (say) of an initial spacelike hypersurface S, which extends to spa- 
tial infinity. Denoting by hij the intrinsic spatial metric on S = 1,2,3), 
Cauchy data would, loosely speaking, consist of hij , (j) and the corresponding 
normal derivatives on 5. By contrast, in quantum theory, one typically asks 
for the amplitude to go from an initial configuration such as {hij,(j))i on an 
initial hypersurface E/ , to a final configuration {hij ,(j))F on a final hypersurface 
S F ■ The problem of finding the quantum amplitude should (naively) be com- 
pletely posed, once one has also specified the (Lorentzian) proper-time interval, 
measured orthogonally between the surfaces E/ and Ej? near spatial infinity. 

Much of the 'non-intuitive' nature of quantum mechanics can be traced to 
the 'boundary- value' nature of such a quantum amplitude [13] , as compared with 
the familiar classical initial- value problem. A crucial mathematical aspect of this 
difference, responsible for a good part of the 'non- intuition', is that the classical 
version of the problem of calculating a quantum amplitude, as posed above, 
would involve solving the classical field equations (typically hyperbolic), subject 
to the given boundary data {hij,(j))i^F on surfaces separated near spatial infinity 
by a Lorentzian time interval T . As is well known, a boundary- value problem 
for a hyperbolic equation is typically not well posed. For typical boundary 
data, a classical solution will not exist [14,15]; or, if it does exist, it will be non- 
unique. The straightforward cure for this ill, due to Fcynman [13], is of course 
to rotate the Lorentzian time-interval T into the complex: T \T\ exp{—id) , 
with < 6* < 7r/2 . 

Consider, for a simple example, a Dirichlet boundary-value problem, but 
where one begins with the 'Riemannian' version, for the Laplace equation (rather 
than the wave equation) in Euclidean 2-space , with data given on an initial 
boundary {r = 0} and on a final boundary {r = T} . Thus, one considers a 



2 



scalar field (j){T ,x), obeying 

^ + ^ = 0; 0<r<T, -oo<a;<+oo, 1.1 

on the assumption that (}> decays rapidly as |.t| ^oo . Suppose, for simplicity, 
that the Dirichlet boundary data arc taken to be of the form 

(j>{T = Q,x) = 0, </)(t = T,x) = M^)- (1-2) 

This problem is treated simply, by taking a Fourier transform with respect to 
X . Define, for example, 

1 r°° 

$i(fc) = ^ e-''''' Mx) dx . (1.3) 



/2tt 

Then the (unique) solution is given by 



1 r .^k. sinh(fcr) 



As with a solution of any elliptic partial-differential equation with analytic co- 
efficients, (j) is automatically a (real- or complex-) analytic function of both 
arguments t and x [14]. 

One may then, as above, rotate the 'time-interval' T into the complex: T 
Texp{i6) , where < ^ < 7r/2 . The integral expression (1.4) continues to give 
the (unique) solution to the differential equation, where it is understood that 
the boundary data (1.2) are left invariant as 9 is varied. This follows since 
(provided < ^ < w/2 ) the denominator in the integrand of Eq.(1.4): 

sinh(A;Te-^^) = cos(fcTsin6') sinh(A;rcos6') - i sin(A;rsin^) cosh(fcrcos6') , 

(1-5) 

is non-zero for all real fe ^ , whence the integrand is smooth for all real k . 

This good behaviour of our linear boundary-value example, in the complex 
case < 9 < tt/2, also follows from more general arguments [16]. For given 
9 € (0,7r/2), let us define a new 'time' coordinate y = Texp{i9), which is 
adapted to the new 'time-interval' Texp{i9) . In terms of the new coordinates 
{y ,x) , the Laplace equation (1.1) reads 

and the boundary conditions read 

cj>{y = Q,x) = 0, cp{y = Te'\x) = M^) • (1-7) 

The potential (t>{y,x) is a complex solution of Eq.(1.6), which is a strongly 
elliptic partial differential equation in the sense of [16] . The property of strong 
ellipticity guarantees existence and uniqueness in this linear example. 
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But, in the Lorcntzian case 9 = 7r/2 , the denominator becomes — i sin(fcT) , 
which has zeros at k = mr/T , for n any integer, positive or negative. Typical 
boundary data 0i(a;) will not have $i(n7r/T) = for a single value of n (n 
integer, n ^ 0). Following this argument, one can show that there is, in general, 
no solution to the Lorentzian (wave-equation) boundary-value problem above. 
Thus, in this simple example, one can already see that the Lorentzian Dirichlet 
boundary-value problem is badly posed. 

In our coupled non-linear gravitational/scalar-field example, the extreme 
case 6 = ^/2 would correspond to a purely Euclidean time- interval |r| , and 
classically one would then be solving the field equations for Ricmannian gravity 
with a scalar field (j) . Since these field equations are 'elliptic modulo gauge' 
- see [16] - one would expect to have a well-posed boundary-value problem, 
with existence and uniqueness. The intermediate case < < 7r/2 requires 
the interval T and any classical solution to involve the complex numbers non- 
trivially. If the problem turns out to be strongly elliptic, up to gauge, then the 
complex case < < 7r/2 would again be expected to have the good existence 
and uniqueness properties of the real elliptic case. 

In practice, one typically treats the case in which the gravitational and scalar 
initial and final data are close to spherical symmetry. Hence, as a leading ap- 
proximation, one begins by studying the spherically-symmetric Einstein-scalar 
system. This is treated in [17] for Lorentzian signature and is outlined in [18] 
for Riemannian signature. In the Ricmannian case, the metric is taken (without 
loss) in the form (3.5,6) below, involving two functions a, b, which depend on 
two coordinates r , r . For Riemannian signature, the field equations are given 
in Eqs.(3.7-ll) below, as partial differential equations for the metric functions 
a, b and the scalar field (p , with respect to r and r . Even in the spherically- 
symmetric case, very little is known rigorously about existence and uniqueness 
for the Riemannian (or complex) boundary-value problem. For this case, nu- 
merical investigation of the weak-field Riemannian boundary- value problem was 
begun in [18], and has recently been extended towards the strong-field region 
[19]. For weak scalar boimdary data, global quantities such as the mass M and 
Euclidean action / appear to scale quadratically, in accordance with analytic 
weak-field estimates [19]. In the limit of strong-field scalar boundary data, it 
may be that a typical pattern will emerge numerically for the general 'shape' of 
the classical Riemannian gravitational and scalar fields. In that case, it might be 
possible to find analytic approximations for the strong-field limit (quite different 
from those valid in the weak- field case), which could provide further analytical 
insight into the solutions of the coupled Riemannian boundary- value problem. 
In particular, it would be extremely valuable to have strong-field approxima- 
tions which were valid into the complex region, with < < 7r/2 . One might 
conjecture that, as one approaches the Lorentzian limit ^ ^ 0+ , for very strong 
spherically-symmetric boundary data, the solutions correspond to classical Ein- 
stein/scalar solutions which form a singularity, surrounded by a black hole. 

Feynman's +ie proposal [13] for computing quantum amplitudes corresponds, 
in our description, to calculating an amplitude (see below) for a complex time- 
separation ]r]exp(— j^), with < ^ < 7r/2, and then taking the limit of the 
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amplitude as ^ 0+ . 

In the case of quantum amplitudes for Lagrangians with Einstein gravity 
coupled to matter, with anisotropic boundary data posed in 'field language', 
such as the case {hij,(f)i_p above, then at least in the asymptotically-flat case 
with time-interval T, one is inevitably led to consider the complex boundary- 
value problem, with T — > |r|exp(— i^) , but with {hij,<p)i^F unchanged. Even 
for fairly small 6 , solution of this boundary- value problem is expected to smooth 
out variations or oscillations of the boundary data, when one moves into the 
interior by a few multiples of the relevant wavelength. If the problem is genuinely 
strongly elliptic, up to gauge, then one will be able to extend the classical 
solution analytically into the complex. 

This and the following papers will use this construction to study, in par- 
ticular, a model concerning nearly-spherical collapse to a black hole, again 
(in the first instance) with Einstein gravity coupled to a real massless scalar 
field (f), except that, for quantum reasons (see Sec. 2), we consider the simplest 
locally-supcrsymmctric model which contains the bosonic Einstein/matter the- 
ory. For this N = 1 supergravity-plus-supermatter model [20], the scalar field 
(p becomes complex, with a massless spin-^ partner, and the graviton acquires 
a spin- 1 gravitino partner. It is assumed that there is a 'background' Rie- 
mannian spherically-symmetric bosonic solution (7^^ , <&) to the Einstein/scalar 
boundary- value problem, where 7^1^ is the background 4-metric and $ the back- 
ground scalar field (taken to be real). For simplicity, one can assume that, near 
the initial surface S/ , the gravitational and scalar fields vary extremely slowly, 
corresponding to diffuse bosonic matter near S/ . For further simplification 
in the quantum calculation, one can assume in the anisotropic case that the 
spatial restriction gij of the exact metric , together with the scalar field 
, is nevertheless spherically symmetric on the initial surface Sj {t = 0) ('no 
incoming particles at early retarded times'), although anisotropic on the final 
surface Ej? . In calculating the Lorentzian quantum amplitude, one would like 
to take at a sufficiently late time T that all the quantum radiation due to the 
evaporation of the black hole will by then have been emitted. In the Rieman- 
nian or complex version, this corresponds to choosing final data {hij , </>) p which 
are nearly spherically symmetric, allowing only for a distribution of weak-field 
(anisotropic) graviton and scalar data on that part of S^? (say, Rq < r < Ri , for 
some large radii Ro ,Ri), which corresponds to the arrival of radiation at Si? . 
The classical back-reaction of the radiation on the geometry can be described 
as follows: on Si?, there is a slowly- varying mass function m{r), where r is 
an intrinsic radial coordinate. The radial rate of change m! (r) is given via the 
Einstein equations, in terms of the averaged energy output in radiation [21,22]. 
Here, m(r) will be extremely small in a region around r = , will then increase 
very slowly through the radiation region, and will then settle at Mj for large 
r , where Mi denotes the conserved ADM (Arnowitt-Deser-Misner) mass of the 
space-time [23,24], as measured, say, at the initial surface S/ . 

In the locally-supersymmetric models of = 1 supergravity coupled to su- 
permatter, such as those in [20,25], the quantum amplitude for a particular 
configuration of bosonic final data, posed on the surface Si?, is expected to 
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be proportional to cxp(— /) or exp(iS) , possibly including finite loop correc- 
tions (see Sec. 2) [15,26,27]. Here / and S are, respectively, the Euclidean and 
Lorentzian classical action corresponding to the solution of the above boundary- 
value problem, with given spherically-symmetric initial data (say) on S/ . Non- 
zero fermionic boundary data and resulting classical actions can also be included, 
provided that the boundary data are suitably posed [28-30] . The resulting clas- 
sical fermionic fields and action will then be elements of a Grassmann algebra, 
as usual in the holomorphic representation for fermions [31]. (In Lorentzian 
signature, this point of view has also been taken in [32] , concerning the Cauchy 
problem for A'^ = 1 supergravity.) Of course, for this statement to make sense, 
one must again write T = \T\ exp{—i9) and take the limit as 6* ^ 0-|- . As 
follows from Sec. 5 and is described in detail for scalar boundary data in Paper 
II, this gives a Gaussian form for the weak-field amplitude. 

The above description refers to the general unitary evolution of states in the 
quantum field theory, assuming that the theory does indeed have meaningful 
quantum amplitudes - this latter point is discussed further in Sec. 2, particu- 
larly in relation to local supersymmetry. Thus, from this viewpoint, there is 
no question of loss of quantum coherence or of information. Prom 1976 until 
July 2004, the most generally accepted option for the end-point of gravitational 
collapse to a black hole was, in fact, that quantum coherence or information 
would be lost [2,6]. The somewhat abrupt change in opinion since then [10] 
now makes it possible at last to begin publishing our late-1990's work on an 
alternative option, outlined above in this Introduction, namely, that there are 
quantum amplitudes (not just probabilities) for final outcomes, and that the 
end- state is a combination of outgoing radiation states. Two letters describing 
the basic outline of this work have appeared or will shortly appear [11,33]. 

The present paper is concerned with setting up some of the underlying lan- 
guage and results that will be needed for further applications (see below). In 
Sec. 2, as mentioned above, we discuss the conditions under which amplitudes in 
quantum field theory (incorporating Einstein gravity) are meaningful: almost 
certainly, this restricts one to theories invariant under local supersymmetry, that 
is, to supergravity or to suitable models of supergravity with supermatter. It 
is expected that such models give finite amplitudes, in the present description 
[15,26,27]. In Sec. 3, we consider the separation of the gravitational and scalar 
fields into spherically-symmetric background parts (7^^,/ , $) as above, and non- 
spherical perturbations at linear and higher orders. The classical field equations 
are considered. In particular, the (classical) back-reaction is considered, in which 
terms up to quadratic in the perturbed metric and scalar field provide an extra 
effective energy-momentum source for the gravitational field. High-frequency 
averaging is also described; this smoothes out the 'random' effects of the emit- 
ted black-hole radiation, to give a simpler treatment of the coupled gravity-plus- 
radiation system. Sec. 4 is concerned with the decomposition of the perturbed 
part of the scalar field into angular harmonics. Following this decomposition, 
the resulting (r , r)-dependent classical scalar field equation is considered; this 
will be essential in the calculation of scalar-field amplitudes in Paper II. In Sec. 5, 
we study the action of classical solutions of the Einstein/massless-scalar field 
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equations; this again will be needed in Paper II. Sec. 6 contains the Conclusion. 

Here, we also describe briefly the main content of this project as presently 
envisaged. Paper II is concerned with the calculation of the quantum amplitude 
for perturbations on the final surface T,p which are purely spin-0 , whereas, for 
simplicity in that calculation, the gravitational field on Si? is taken to be ex- 
actly spherically symmetric. Next, we will relate the description of Paper II to 
the Bogoliubov transformations [33] which arc familiar from many earlier treat- 
ments of particle creation [3,34,35]. In the course of this, the relation between 
the Vaidya metric [21] and the space-time geometry in the region containing the 
outgoing radiation flux will be worked out carefully. The quantum amplitudes 
or wave functionals for our 'complex-rotated-T' problem, as in Paper II, are 
further related to an alternative description in terms of coherent and squeezed 
states; this points to a more universal validity of the procedure adopted in this 
project, with further applications to quantum amplitudes for inhomogeneities 
in cosmology, for example. Finally, the procedure of Paper II for spin-0 per- 
turbations will be generalised to the case of spin-1 Maxwell (or Yang- Mills) 
perturbations, and to spin-2 graviton perturbations. The fermionic case s = 5 
of massless neutrino perturbations will be treated similarly. The remaining 
fermionic case of the spin-| field, needed for locally supersymmetric models, is in 
preparation. In addition, substantial work has been done concerning computer 
solution of the elliptic or the complexified version of the spherically-symmetric 
Einstein/massless-scalar field equations, following [18]. Results of this further 
work should soon be presented [19]. 

2 The quantum amplitude for bosonic boundary 
data 

Consider the 'Euclidean' quantum amplitude to go between prescribed initial 
and final purely bosonic data, each given on a 3-surface which is 'topologically' 
(diffeomorphically) M.^ , and each carrying an asymptotically-flat 3-mctric. It 
is further necessary to specify the proper (Euclidean) distance r , measured or- 
thogonally between the two surfaces at spatial infinity. This amplitude will be 
given (formally, at least) by a Feynman path integral over all Riemannian infill- 
ing 4-geometries together with any other bosonic fields, each such configuration 
being weighted by exp(— /) , where I is the 'Euclidean action' of the configura- 
tion. If this definition is meaningful, one expects that the resulting 'Euclidean' 
quantum amplitude has the semi-classical form 

Amp ~ {Ao + hAi+h'^A2 + ...) exp{-lB/h) , (2.1) 

asymptotically in the limit that [Ib /ti) . Here. Is is the classical 'Euclidean 
action' of a Riemannian solution of the coupled Einstein and bosonic-matter 
classical field equations, subject to the boundary conditions. For simplicity, we 
assume that there is a unique classical solution, up to gauge and coordinate 
transformations. But it is quite feasible, in certain theories and for certain 
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boundary data, to have instead (say) a complex- conjugate pair of classical so- 
lutions [36]. 

The classical action Ib and loop terms Aq ,Ai ,A2 . . depend in principle 
on the boundary data. In the case of matter coupled to Einstein gravity, each of 
Ib , Ao , Ai , . . . will also obey differential constraints connected with the local 
coordinate invariance of the theory, and with any other local invariances such 
as gauge invariance (if appropriate) [15,37]. 

In particular, when the theory is also invariant under local supersymmetry, 
the semi-classical expansion (2.1) may become extremely simple [15,26,27]. For 
example, fov N = 1 supergravity, one has [15]: 

Amp ~ Ao e^p{-I/h) . (2.2) 

In this theory, the one-loop factor is in fact a constant. When one allows the 
boundary data to include both bosonic and fcrmionic parts, suitably posed, one 
expects that a classical solution of the coupled bosonic/fermionic field equations 
will still exist. In this case, I in Eq.(2.2) denotes the classical action, includ- 
ing now both bosonic; and fcrmionic contributions. Related properties hold for 
= 1 supergravity coupled to gauge-invariant supermatter [25,26]. There will 
also, for example, be analogous consequences for locally-supersymmetric mod- 
els oi N ^ 1 supergravity coupled to supermatter [26,27], which in the simplest 
case include a complex scalar field and a spin-i field [20]; one expects the ampli- 
tude semi-classically to be of the form (2.1), with finite loop terms Aq ,Ai , . . . 
[15,26,27]. 

In the case (2.2) of = 1 supergravity, the classical action is all that is 
needed for the quantum computation. A corresponding situation arises with 
ultra- high-energy collisions, whether between black holes [38], in particle scat- 
tering [39], or in string theory [40]. 

In the asymptotically- flat, spatially-M^ context appropriate here, the purely 
Riemannian case above corresponds, in Lorentzian-time language, to a (Lorcntzian) 
time-separation at spatial infinity of the usual rotated form T = —ir, where r is 
the (positive) imaginary-time separation defined above. If the four-dimensional 
classical bosonic part of the solution is to be real, then certainly the bosonic 
boundary data should be chosen to be real. Following the standard route, one 
should then study the (now complex) bosonic amplitude (2.1) or (2.2), as T 
is rotated through a range of angles 6 , starting from 6 = 7^/2 and ending at 
e =+e (> 0) , with 

T = T exp{-ie) . (2.3) 

Provided that there continues to exist a (complex) bosonic classical solution to 
the Dirichlet problem, as one rotates 9 from tt/2 towards zero, the expression 
(2.1) or (2.2) should continue to give the form of the quantum amplitude, which 
should be analytic in 6 (among other variables). In particular, this would oc- 
cur if strong ellipticity [16] held for the coupled Einstein/bosonic-matter field 
equations, up to gauge. 
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3 The approximate 4-dimensional metric 



The classical background bosonic fields are at present taken, for simplicity, to be 
only the metric and the masslcss scalar field (f> . In later work, we shall study 
cases in which other-spin fields, with s = ^ , 1 or | , are included as perturba- 
tions of a spherically-symmetric background solution [30,41,42]. The classical 
solutions (g^i, , (p) of the coupled Einstein/scalar field equations below are taken 
to have a 'background' time-dependent spherically-symmetric part {"ffn, , <&) , 
together with a 'small' perturbative part (/i^i/ , (/'pert)- The perturbative fields 
hfii, and 4'pert , which live on the spherically-symmetric background 4-geometry 
j^i, , can, as usual, be expanded out in terms of sums over tensor (spin-2), vector 
(spin-1) and scalar harmonics [43,44]. Each harmonic is weighted by a function 
of the Riemannian time- and radial coordinates (t, r) . 
The Einstein field equations are taken in the form 

Gij,v = Rjxv — -Rg^v = SttT^k , (3.1) 

where R/jtv denotes the Ricci tensor, R the Ricci scalar and T^^ the energy- 
momentum tensor. For a (real) massless scalar field (j) , one has 

T^u = </>,M</',.-^5M-^(0,a0,/3<7"O . (3.2) 

The gravitational field equations further imply the scalar field equation (the 
Laplace-Beltrami equation [45] ) : 

9^.(9'/^ 9^" = , (3.3) 

where g denotes det(5(^,^) , and at present we assume that g^i, is real Rieman- 
nian, whence <? > . 

The corresponding Riemannian variational principle refers to an action func- 
tional of the form [15] 



/ = ^ 
167r 



J d*x g^/"^ R + ^ y "^^2; g^^'^{y4>f + boundary contributions . 

(3.4) 

The appropriate boundary terms will be discussed below in Sec. 5. 

In the Riemannian case [18], the 'large' or 'background' 4-metric can be put 
in the form: 

rfs^ = e'' dr^ + e^dr"^ + r^{de^ + sin^O d<tP-) , (3.5) 

where 

b = b{T,r) , a = a{T,r) . (3-6) 

If the gravitational field were exactly spherically symmetric, as in Eq.(3.5), and 
if the scalar field were also rotationally invariant, being of the form 0(r, r), then 
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the Ricmannian sphcricaUy-symmctric scalar and Einstein field equations would 
hold [18]. The scalar field equation reads: 

^ + e^-" cP" + ^{d-b) ^ + r-^ e''-" (l + e") = , (3.7) 

where (') denotes d{ )/5t and ( )' denotes d{ )/dr . Together with Eq.(3.7), a 
slightly redundant set of gravitational field equations is given by: 

a' = - 47r r (e""'' - (/)'^ ) + r"! (l - e"") , (3.8) 
b' = - 47r r (e"-'' _ ^'^^ _ (l - e") , (3.9) 
a = 8Trr<j)(j)', (3.10) 

a + e^-'^b" + ^{a-b)d - r-^e''-''{l-e''){b' + 2r-^) = Stt (^^ ^ e^"" </.'') . 

(3.11) 

The metric and the classical field equations in Lorentzian signature [17], or 
for certain types of complex metrics, can be derived from the above by the 
formal replacement 

t = r e-^^ , (3.12) 

where 9 is independent of 4-dimensional position, and where should be rotated 
from to 7r/2 . 

In the bosonic black-hole evaporation problem, the classical Riemannian 
metric and scalar field will not be exactly spherically symmetric; similarly for 
any non-zero spin-i and spin-| classical (odd Grassmann-algebra- valued [15]) 
fermionic solutions in the locally-supersymmetric generalisation [20] . In particle 
language, rather than the field language which is mostly being used in this paper, 
huge numbers of gravitons and scalar particles will continually be given off by 
the black hole (together with any fermions allowed by the model), leading to 
an effectively stochastic distribution, in which, for any given spin s , the field 
fluctuates around a spherically-symmetric reference field. 

Consider, for example, gravitational and scalar perturbations about a Rie- 
mannian spherically-symmetric reference 4-metric 7^jy and reference scalar field 
$ . In perturbation theory in general relativity [46] , one considers a one- 
parameter (or many-parameter) family of 4-metrics, here of the form 

5^,(x,e) = 7^,(x) + e/iW(a;) + /i(fj(x) + ... , (3.13) 

where h^}} is the first-order metric perturbation, /i^^J is the second-order per- 
turbation, etc. Throughout, the superscript will refer to the background, 
while denotes terms linear in perturbations, etc. Indices are to be raised 
and lowered using the background metric 7^'',7^i/ . Covariant derivatives with 
respect to the background geometry are denoted either by a semi-colon -a or 
equivalently by Va • 

Analogously, we split a real massless scalar field (j) into a spherically-symmetric 
background piece 3>(t, r) and a (non-spherical) perturbation: 

(t>{x,e) = $(r,r) + e <f)^^\x) + e"^ (t>^'^\x) + ... . (3.14) 
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The spherically-symmetric background part $ will be non-zero if the back- 
ground scalar data (j) at early and late Euclidean times r contain a non-trivial 
spherically-symmetric component. The perturbation fields ^'^^(a;), (j)^'^\x), ... 
will, in general, contain all non-sphcrical angular harmonics. These fields must 
be chosen such that the entire coupled Einstein/scalar system satisfies the clas- 
sical field equations, as well as agreeing with the prescribed small non-spherical 
perturbations in the initial and final data, both gravitational and scalar. The 
effective energy-momentum source for the spherically-symmetric part of the 
metric includes contributions formed quadratically from the non-spherical grav- 
itational and scalar parts and — see Eqs. (3.27-29) below for further 
detail. 

In the simplest case, one can restrict attention to the exactly spherically- 
symmetric Ricmannian model of Eqs. (3. 5-11). Then the background metric 7^jy 
will correspond to the metric (3.5,6), with respect to suitable coordinates, and 
#(t, r) here will correspond to the 0(r, r) of those equations. This Riemannian 
boundary- value problem, involving a system of coupled partial differential equa- 
tions in two variables (r, r), has been studied numerically in [18] for particular 
choices of boundary data, and is currently being investigated in much greater 
detail [19]. 

By contrast, the Lorentzian-signature version of the spherically-symmetric 
classical Einstein/scalar system must be studied as an initial- value evolution 
problem, in order to be well-posed [17]. One conceivable initial profile for the 
scalar field, which has been much studied in the Lorentzian-signature numerical 
problem [47,48], is an ingoing 'Gaussian' shell of scalar radiation. To define such 
initial data, work in a nearly-flat space-time at very early times (large negative 
Lorentzian time-coordinate t). Define an advanced null coordinate 



where fo ,k ,d,ro and A are all positive real parameters. Here, the radial 
extent, Lq , of the 'Gaussian' is given by Lq ^ A . The numerical evolution of 
such initial data provides a model of spherical collapse. In particular, two main 
qualitatively different regimes of initial data can be distinguished. Firstly, if Lq 
or A is too large, then the initial data are 'diffuse', there is little self-interaction, 
and the incoming scalar profiles pass more or less straight through each other, 
leaving behind nearly-flat space-time plus small perturbations. Secondly, if A 
(or Lq) is less than a certain critical value, the interaction is sufficiently non- 
linear that a black hole forms. 

Returning to the Riemannian or to the complex case, one can expand out 
the Einstein field equations (3.1,2) in powers of e . At lowest order (e°), one has 



V = t + r 



(3.15) 



The incoming 



'Gaussian' shell is asymptotically, at early times. 




(3.16) 
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the background Einstein and scalar field equations 

<^ - \r^'^1,^ = SttTW , (3.17) 
7''"^;^. = , (3.18) 

where R^^} denotes the Ricci tensor and i?^''^ denotes the Ricci scalar of the 
background geometry 7^^ , and where a semi-colon now denotes covariant dif- 
ferentiation with respect to the background geometry. Further, 

tW = _ l^^^($,^$,^^a/3) (3.19) 

denotes the background spherically-symmetric energy-momentum tensor. These 
field equations are equivalent to Eqs.(3.7-ll), when coordinates are taken as in 
Eqs.(3.5,6). 

The linearised (e^) part of the Einstein equations reads (see Section 35.13 of 
[24]) 



Here, h'j},^ is defined by 



(3.20) 



h'-^} = h^^^-l^^uh^'^ , (3.21) 

where 

Here, R^y^ii^a denotes the Riemann tensor of the background geometry 7^^!/ . 
Further, Tj^D denotes the Unearisation or O(e^) part of the energy-momentum 
tensor T^i^(x, e). Explicitly, 

tW = 2V(^<^« V.)$-7^.(V„*)(V"(/.«) + ^ (7^./i«"''-ft(,l)7"'')(V.$)(V, 

(3.23) 

The linearised Einstein equations (3.20-23) are most easily studied in a 'lin- 
earised harmonic gauge' [24,49] , in which, by an infinitesimal coordinate trans- 
formation, one arranges that 

h'l^l''" = (3.24) 

everywhere. Since the gravitational background 7^^ is spherically symmetric, 
the linearised Einstein equations (3.20-23) can be further decomposed into three 
independent sets of equations. These describe repectively scalar (spin-0) per- 
turbations associated with matter-density changes (T^r'), vector (spin-1) per- 
turbations associated with matter- velocity changes (T^l^), and gravitational 
radiation (spin-2) associated with anisotropic stresses [T^p) [50]. 
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The linearised (e^) part of the scalar field equation (3.3) yields 

j^'-cPlll -{U'^^''^,,),^ = . (3.25) 

Thus, the hncariscd Einstein and linearised scalar field equations (3.20-23,25) 
are coupled. 

At order , the gravitational field equations give the second-order contri- 

(2) 

bution Gil,; to the Einstein tensor 

G^iv = Rfjtv — -^Rg/xu ; (3.26) 

details are given below. We note that Gl^J includes a well-known contribu- 
tion quadratic in the first-order perturbations h'j}} and their derivatives — see 

Eq. (35.58b) of [24]. This part represents an effective energy-momentum-stress 

(2) 

density due to the gravitational perturbations, including gravitons. Gfj.^ also 
contains contributions at quadratic order, formed from the background $ and 
the linearised (p^^^ together with 'Jij^i, and /ij/j . These parts represent the O(e^) 
contribution of the scalar- field energy- momentum tensor T^^ of Eq.(3.2). 

Explicitly, one finds, after a lengthy calculation [9], that the Einstein equa- 
tions, to quadratic order in perturbations, read 

= 8^T^°) + 8TT T^'J + SttT;, - . (3.27) 

Here, 

(3.28) 

and 

I (OCT 

(3.29) 

where the explicit form of will not be used here. The above expressions 
will be needed particularly in studying the Vaidya metric [21], which describes 
approximately the late-time region of the geometry following gravitational col- 
lapse to a black hole, containing a nearly-steady outgoing flux of radiation. The 
Einstein field equations, averaged over small regions, give the contribution of 
massless scalar particles, gravitons, etc., to the nearly-isotropic fiux. 
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Physically, for the Ricmannian or complex boundary-value problem dis- 
cussed in Secs.l and 2, the 0(e) perturbations in the 4-metric g^i, and scalar field 
0, relative to the spherically-symmetric background solution (7^^ , $), should 
arise classically from 0{e) perturbations away from spherical symmetry in the 
boundary data gij , $ (or d^/dn) at the initial and final surface. Moreover, as 
mentioned earlier, provided that the perturbed boundary data contain numer- 
ous high harmonics, the 4-dimcnsional perturbations in the interior would be ex- 
pected to have an effectively stochastic nature. When averaged over a number of 
wavelengths, the effective perturbative energy-momentum tensor above, T^J^^, 
will yield a spherically-symmetric smoothcd-out quantity < T^J^^ > [51,52]. 
In the locally-supersymmetric version of this theory, the energy-momentum ten- 
sor due to the spin-i and spin-| fields will also contribute to < T^f^ > . This 
will account, in particular, for the gradual loss of mass by radiation of a black 
hole in the nearly-Lorentzian sector (that is, in the case of a time-interval at 
infinity of the form T = r exp(— i^) , where r is real and positive, and ^ = e is 
small and positive). Although < T^f^ > is small, being of order e^, its effects 
on the black-hole geometry, including the mass, will build up in a secular fash- 
ion, over a time-scale of order 0(e~^). Such secular behaviour appears often in 
perturbation problems [53,54] — for example, in the familiar treatment of per- 
ihelion precession for nearly-circular orbits in the Schwarzschild geometry [55]. 
In our boundary- value problem, whether regarded as classical or quantum, the 
initial boundary data will be spread over a 'background' extent of 0(1) in the 
radial coordinate r on the initial surface T,j . But, corresponding to the 0(e~^) 
time-scale for the black hole to radiate, the final data on will be spread 
over a radial-coordinate scale of O(e^). Thus, even the classical boundary- value 
problem here is an example of singular perturbation theory [53,54]. 

The standard treatment of high-frequency averaging in general relativity was 
given by Brill and Hartlc [51] and by Isaacson [52]. Let < > denote an average 
over a time Tq much longer than typical wave periods, together with a spatial 
average over several wavelengths A . Then 

< > = 7/.-^ , < <^ > = $ , (3.30) 

<</)(^)> = , </iW> = , (3.31) 

<a<,/iW> = , <d„dph<;}}> = . (3.32) 

Indeed, 

< C(0) > = , (3.33) 

for any background quantity (7^°^ . Rules for manipulating these averages in the 
high-frequency aproximation are set out in [52]. Under integrals, the average of 
total divergences can be neglected. For example, 

< h^''"-, h^l;.. > = - < h^'\ C > • (3-34) 

Further, covariant derivatives commute for high-frequency waves. The rules 
(3.30-34) imply 

<r«> =0 , (3.35) 
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<T^2;> = <V(^</.«V.).^«-^7/..V„<^«V"</)W > . (3.36) 

We can now rewrite the background field equations (3.7-11) in a form smoothed 

out by averaging over a number of wavelengths of the scalar and gravitational 
perturbations [51,52]. The equation which includes the quadratic-order contri- 
bution of the perturbations as a source for the background geometry reads: 

G(°j (7) = 8n T(°) + Stt (< T^l^ > + <T'^^>) . (3.37) 

The terms in this equation vary over length-scales >> A . The 'source equation' 
for h^^v , analogous to Eq.(3.20) for h^}} , is 

G«(7,/^('^) = 87r(T;,-<T;,>) + 8^(T(2)-<T(2)>) . (3.38) 

Here, the left-hand side G^V^ (7 , h^"^^ ) denotes the first perturbation of the Ein- 
stein tensor Gfj,^ about the background metric , but with its linear argument 
taken to be h'^} rather than . Thus, — 2G|/j(7,/i(^^) is given by the left- 
hand side of Eq.(3.20), subject to Eqs.(3.21.22). Hence, the left-hand side of 
(2) 

Eq.(3.38) is linear in hJiJ and its derivatives, whereas the right-hand side is 
quadratic in first-order fluctuations. By contrast with Eq.(3.37), the terms in 
Eq.(3.38) vary over length-scales of order A . 

4 Scalar field: harmonic description 

Consider small bosonic perturbations (t>^^'' and h^^} , obeying the linearised clas- 
sical field equations (3.20) and (3.25) about a spherically-symmetric classical 
solution (<i> , j^u) of the Riemannian field equations (3.7-11) for Einstein grav- 
ity coupled minimally to a massless scalar field. The background spherically- 
symmetric data for $ and 7^1/ are posed, as in Sees. 1,2, on the initial and final 
3-dimensional boundaries, separated at spatial infinity by a 'Euclidean time- 
separation' r > . Similarly, the linearised classical perturbations ^^^^ and 
are to be regarded as the solutions to a coupled linear elliptic problem, subject 
to prescribed linearised perturbations (p^^^ (say) and h^j^ on the initial and final 
boundaries. 

Because of the spherical symmetry of the background ($,7^^), one may 
expand the Riemannian 4-dimensional perturbation (/f*-^' in the form 

^ 00 i 

0W(r,r,^,<^) = - ^ ^ F,„(n)i?,„(r,r) . (4.1) 

e=0 m=-i 

Here, Ygmi^) denotes the {i,m) scalar spherical harmonic of [56]. 

Similarly, a generic Riemannian metric pertubation h^/}} may be expanded 
out as a sum over spin-2 (tensor), spin-1 (vector) and spin-0 (scalar) {£,m) 
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harmonics, each weighted by a function of r and r [43,44,57-60]. The amph- 
tudes for graviton (spin-2) emission following black-hole collapse will be treated 
subsequently, including further details of the spin-2 harmonics. But note that, 
because of the coupled nature of the linearised field equations (3.20,25) for cj)^^^ 
and /i^i^ , the resulting linear field equations in r and r for Rim{T, r) of Eq.(4.1) 
and its gravitational analogues will also be coupled in the strong-field 'collapse' 
region of the 'space-time'. 

The boundary conditions on the radial functions R£m{T ,r) as r ^ follow 
from the regularity there of the whole Riemannian solution, consisting of and 
the 4-mctric g^j^ (but viewed in 'ncarly-Cartcsian coordinates' near r = 0). This 
regularity of the solution in turn follows since the coupled field equations are 
'elliptic modulo gauge'. For simplicity, the boundary data, on both the initial 
and the final 3-surfac;c;, must be chosen to be sufficiently regular or smooth over 

, in addition to being asymptotically flat. Even when one takes a complex 
Lorentzian time-separation-at-infinity 



as in Eq.(2.3), with < 9 < 11/2 , one expects that the field equations (up to 
gauge) will be strongly elliptic [16], whence all classical fields must be analytic 
in the interior of the large cylindrical boundary formed by the initial and final 
surfaces, together with a surface at large r . 

Suppose that the boundary conditions on the final surface are taken to con- 
sist of weak and very diffuse scalar and gravitational fields, regarded as pertur- 
bations of fiat 3-space . (One also requires that the ADM (Arnowitt-Deser- 
Misner) mass of the final intrinsic boundary 3-metric gij , computed from the 
(1/r) part of the fall-off of gij to the flat metric Sij [23,24,61] should be the same 
as the ADM mass of the initial surface. This will be discussed further in Sec. 5 
below.) Physically, such weak and diffuse final boundary data may be imagined 
to be a possible late-time remnant of gravitational collapse, namely a snap-shot 
of a large number of scalar particles and gravitons, as they make their way out 
to infinity. Near the final surface, then, the coupling in Eqs. (3.20,25) between 
the linearised perturbations (j)^^^ and will almost have disappeared. The 
perturbed scalar field equation at late times is simply 



with respect to the spherically-symmetric background geometry 7^,^ . 

Making the mode decomposition (4.1) of (f>^^\ one obtains the {i,m) mode 
equation 



T = T exp{—i9) 



(4.2) 








(4.3) 



+ ^ {dr{a - b)) {drRem) ' Vi{T , r) Ri, 



tm 



= 



(4.4) 



Here, 



Ve{T,r) = 



£{£ + 1) + 



2m (r , r) 



(4.5) 



r 
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and m{T ,r) is defined by 

exp(-a(T,r)) = i _ Mill) . (4.5) 

In an exact Schwarzschild solution with no scalar field, one would have e*" = 
e~'^ = 1 — ^ , with M the Schwarzschild mass; in that case, m(r, r) would be 
identically M . The potential V^(T,r) of Eq.(4.5) generalises the well-known 
massless-scalar effective potential in the exact Schwarzschild geometry [24], 
which vanishes at the event horizon {r = 2M} and at spatial infinity, and 
has a peak near {r = 3M}. 

The definition (4.6) of m(T ,r) is also consistent with the usual description 
of the Lorentzian-signature Vaidya metric [21,22]. In terms of a null coordinate 
u and an intrinsic radial coordinate r , the Vaidya metric reads 

ds^ = - 2 du dr-(^l- ^HM^ du" + [dO^ + sin^ 9 dcj)^) . (4.7) 

Here, m{u) is a monotonic-decreasing, but otherwise freely specifiable smooth 
function of ?i, corresponding to a suitable spherically-symmetric outflow of null 
particles, for example the energy-momentum tensor of a black hole evaporating 
via emission of scalar particles at the speed of light. The Vaidya metric has 
been used often to give an approximate gravitational background for black- 
hole evaporation at late times [4,5,62]. This connection will be treated more 
thoroughly in our subsequent work. 

There is, of course, an analogous decoupled harmonic decomposition, valid 
near the final surface, for the weak gravitational-wave perturbations about the 
spherically-symmetric background - again described in subsequent work. For 
simplicity of exposition, we shall in the following Paper II restrict attention to 
weak- field final configurations for spin-0 (scalar), and calculate their quantum 
amplitudes on the further assumption that the final 3-metric hijp is exactly 
spherically symmetric (in addition to the assumed spherical symmetry of the 
initial data (f)i and hiji). Once the methods are established in the simplest 
spin-0 case, generalisation to the case of higher-spin fields becomes relatively 
straightforward. 



5 The classical action 

Consider, for definiteness, an asymptotically-flat Lorentzian-signature classical 
solution {gij,v,(j>) of the coupled Einstein/massless-scalar field equations, be- 
tween an initial hypersurface S/ and a final hypersurface Si? , separated by a 
proper Lorentzian time T at spatial infinity. Write S for the Lorentzian action 
functional, which corresponds to the Riemannian action functional / of Eq.(3.4) 
with suitable boundary contributions [15], appropriate to fixing the boundary 
data {hij , <p)i and {hij ,(f>)F , according to iS = —I . Then, at the Lorentzian- 
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signature solution above, one has [15,63] the classical action 



'S'cla 



3277 VJe 



Si 



2 „ 

(5-1) 

Here, tt*-' = tt-'* is (167r) times the Lorentzian momentum conjugate to the 
'coordinate' variable hij on a space-like hypersurface, in a 3 + 1 Hamiltonian 

decomposition of the Einstein/masslcss-scalar theory [64]. Explicitly, in terms 
of the Lorcntzian-signature second fundamental form Kij = Kji of the hyper- 
surface [15,49], TT*-' is given by 

TT^J = h^/^{K^i -Kh'^) , (5.2) 

where h = det(/ijj) and K = h^^ . Further, tt^ is the Lorentzian momentum 
conjugate to the 'coordinate' variable (j) ■ Explicitly, 

TT^ = /i^n^(V^<^) , (5.3) 

where n'^ denotes the (Lorentzian-signature) future- directed unit time-like vec- 
tor normal to the hypersurface. 

Suppose instead that one has a complex or a Riemannian solution (g^^ , cf)) 
between asymptotically-fiat boundary data {hij , </>)/ and {hij ,(j))F on initial and 
final hypersurfaces S/ , S^^- , where the time-separation T at infinity has the form 
T = T exp{—id) , as in Eq.(2.3), where r is positive real and < 9 < Tr/2 . This, 
as above, is expected to provide the most natural arena for asymptotically-flat 
boundary-value problems involving gravitation, if strong ellipticity holds, up 
to gauge. For such a solution, the Lorentzian-signature classical action Sdass 
continues to be defined by Eq.(5.1). This will in general be complex, although 
for a real Riemannian solution with 9 = 7r/2 , the Riemannian action /dass j 
defined by /dass = ^ i Sdass , will be real. The boundary contribution at spatial 
infinity to the Riemannian action functional / corresponding to Eq.(5.1) is +Mt 
[15]. The boundary contributions to the functional /, due to the presence of 
the boundaries S/ and Hp with data {hij , and {hij , (^)f specified on them, 
are 



Ii + If = if/ - / )d'x eTT^^' hij + U [ - I ) d'x eTT^ 



327r 

Here, 



(5.4) 



eTT'^' = ( ^j^v _ h'^) (5.5) 

is given by the same formula as tt'-' in Eq.(5.2), except that Kij has been re- 
placed by the 'Euclidean' second fundamental form e-^ij ) as defined and used 
in Eqs.(2.6.23,24) of [15]. In particular, 

^Kij = - iKij . (5.6) 
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Similarly, the scalar-momentum variable tt^ of Eq.(5.3) has been replaced by its 
'Euclidean' version e7r^ , defined by 

eTT^ = h^'^ ^n^ (V^(/.) , (5.7) 

where [15] 

^ni^ = -in'' (5.8) 
denotes the unit future-directed Riemannian normal. 

The quantity M in Eq.(5.1) denotes the Arnowitt-Deser-Misner (ADM) mass 
of the 'space-time', as measured near spatial infinity from the (1/r) part of the 
fall-off of the intrinsic spatial metric hij on S/ and [23,24]. As mentioned 
in Section 5, it is essential, for a well-posed asymptotically- flat boundary- value 
problem, that the intrinsic metrics hiji and hijp be chosen to have the same 
value of M . Otherwise, if M/ ^ Mp , then any classical infilling 'space-time' will 
have S/ and T^p badly embedded near spatial infinity, and the entire 4-metric 
gfj_,y will not fall off as rapidly as it should, as r ^ oo [61]. 

In applications to black-hole particle emission, we naturally make use of the 
pcrturbative splitting g^^ = 7^,^ + h^/}} + • • • , (j) = ^ + 4'^^'^ + . . . , where the 
spherically-symmetric 'background' (7^1^ , $) obeys the coupled Einstein/massless- 
scalar classical field equations, as does the full classical solution (p^^, , <p) . (The 
formal device of including a small parameter e has been relaxed here: we now 
set e = 1 .) 

The linearised fields /ij/j and (p^^^ may be decomposed into sums of appropri- 
ate angular harmonics labelled by quantum numbers {£ , m) as in Section 5, and 
without loss of generality it may be assumed that any spherically-symmetric 
1 = linear-order perturbation modes have been absorbed into the spherically- 
symmetric background (7^^ , $) . Then (say) the Lorentzian classical action 
'S'ciass of Eq.(5.1) may be split as 

•Sclass = -^ciiss + "^class + '^'class + • • ■ ■ (5-9) 

Here, S'^J^ss ^^c background action, given by Eq.(5.1), but evaluated for the 
spherically-symmetric solution (7^1, , $) . The mass M appearing in S^^^^ will 

(2) 

be that determined from (7ij)/ or {jij)F ■ The next term is S^^^^ , formed 
quadratically from the linear-order perturbations; one may verify that the linear- 
order term S^^^^^ is zero, because of the above definitions. In an obvious notation, 
one has 

= 3^ (i -/J "S" - 5 il -0 ^'^ ^ 

(5.10) 

(2) 

Note that there is no contribution to the second-order expression •S'^jg^gg from the 
—MT term in Eq.(5.1), again because of the above definitions. 

The expression (5.1) for 5'ciass[(^ij ? 0)/ ; (^ij ? 4')f ? together with the asymp- 
totic series (5.9) for the classical action and the expression (5.10) for S^^^^^ 
formed from the linearised perturbations, will be basic in calculations concern- 
ing quantum amplitudes in subsequent work. 
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6 Conclusion 



This paper has been concerned with setting up a basic framework and formaUsm 
for treating quantum amphtudcs involving possibly strong gravitational fields, 
governed by a Lagrangian containing a locally supcrsymnictric version of Ein- 
stein gravity and matter. This includes the case of gravitational collapse to 
a black hole, but is considerably more general, being applicable also to many 
cosmological problems involving small fluctuations of an isotropic homogeneous 
universe, in the quantum context. 

The underlying approach in this paper has been to calculate the quantum 
amplitude to go from data (both for gravity and any matter fields) specified 
on an initial spacelike hypersurface S/ to corresponding data given on a final 
hypersurface "Ep . Since, in the black- hole context mainly studied here, the 
space-time should be asymptotically flat, we take the simplest case in which 
both E/ and T,p are diffeomorphic to Euclidean space R^, and such that their 
intrinsic 3-dimensional metrics hiji , hijp are asymptotically flat at spatial in- 
flnity. The 'boundary data', which should determine the quantum amplitude 
uniquely, then consist of hij and suitable components of any matter flelds, on S/ 
and E p , together with the (Lorentzian) proper time- interval T between E/ and 
Ei? , measured near spatial infinity. Following Feynman's +i€ prescription [13], 
we rotate T into the complex: T \T\ cxp{—i9) , for < 9 < 7r/2 . One expects 
that the corresponding complex classical boundary- value problem for the Ein- 
stein/matter field equations should be well-posed, unlike the purely Lorentzian 
case with T real. The remaining analysis is mainly concerned with proper- 
ties and consequences of such complex solutions of the field equations, where 
the data h^j , etc., for gravity and matter are held fixed on the boundaries E/ 
and T,p , but with T — > \T\ exp(— i6) . The quantum amplitude for linearised 
perturbations is given principally through the second-variation classical action 

(2) 

■^ciass ' ^ functional of the boundary data, with the amplitude proportional 

(2) 

to exp(i -S'^iass) , except near Planckian energies. Feynman's prescription then 
requires that we take the limit 6^0^ to obtain the Lorentzian amplitude. 

In the following Paper II of this series, we shall evaluate the above ampli- 
tude for a model with Einstein gravity and minimally- coupled scalar field (p , 
in the case that the perturbations on S^- are only in (j) , but not in the final 
gravitational data hij , which are there taken to be spherically symmetric. In 
particular, this is applicable to the case of spin-0 (scalar) radiation from gravi- 
tational collapse to a black hole. In further work, this description will be related 
to the alternative language of Bogoliubov transformations, in which much of the 
earlier work on black-hole evaporation was cast [3,34,35]. Next, we will anal- 
yse in some detail the semi-classical description of the region of the space-time 
containing the outgoing flux of (here, spin-0 and spin-2) radiation, with the 
help of the Vaidya metric [21]. Yet another kind of 'transformation' can also be 
analysed, namely, that between the language of this work and the language of 
coherent and squeezed states; this approach makes it easier to take an overview 
and see (for example) the similarity between our 'local-collapse' or 'black-hole' 
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description and a cosmological version, in which the only boundary is (say) a 
compact 3-surface, such as a 3-sphere [15,36]. The calculations of Paper II 
for the quantum amplitude for purely scalar-field (spin-0) emission will be gen- 
eralised to the cases of spin-1 Maxwell (or Yang-Mills) radiation and to spin-2 
graviton (gravitational-wave) emission. The fermionic case of a massless s = ^ 
(neutrino) field will also be treated. Work on further aspects of the approach, 
as outlined in the Introduction and in [11,33], has also been carried out and 
should soon be completed. 
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